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Abstract 

We investigate analytically and numerically nonlinear vortex spin torque oscillator dynamics in a circular magnetic 
nanodot induced by a spin-polarized current perpendicular to the dot plane. We use a generalized nonlinear Thiele 
equation including spin-torque term by Slonczewski for describing the nanosize vortex core transient and steady 
orbit motions and analyze nonlinear contributions to all forces in this equation. Blue shift of the nano-oscillator 
frequency increasing the current is explained by a combination of the exchange, magnetostatic, and Zeeman 
energy contributions to the frequency nonlinear coefficient. Applicability and limitations of the standard nonlinear 
nano-oscillator model are discussed. 
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Background 

Spin torque microwave nano-oscillators (STNO) are in- 
tensively studied nowadays. STNO is a nanosize device 
consisting of several layers of ferromagnetic materials 
separated by nonmagnetic layers. A dc current passes 
through one ferromagnetic layer (reference layer) and 
thus being polarized. Then, it enters to an active mag- 
netic layer (so-called free layer) and interacts with the 
magnetization causing its high-frequency oscillations 
due to the spin angular momentum transfer. These oscil- 
lation frequencies can be tuned by changing the applied 
dc current and external magnetic field [1-3] that makes 
STNO being promising candidates for spin transfer 
magnetic random access memory and frequency-tunable 
nanoscale microwave generators with extremely narrow 
linewidth [4]. The magnetization in the free layer can form 
a vortex configuration that possesses a periodical circular 
motion driven by spin transfer torque [1,5-11]. For 
practical applications of such nanoscale devices, some 
challenges have to be overcome, e.g., enhancing the STNO 
output power. So, from a fundamental point of view as 
well as for practical applications, the physics of STNO 
magnetization dynamics has to be well understood. 
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In the present paper, we focus on the magnetic vortex 
dynamics in a thin circular nanodot representing a free 
layer of nanopillar (see inset of Figure 1). Circular nano- 
dots made of soft magnetic material have a vortex state 
of magnetization as the ground state for certain dot radii 
R and thickness L. The vortex state is characterized by 
in-plane curling magnetization and a nanosize vortex core 
with out-of-plane magnetization. Since the vortex state 
of magnetization was discovered as the ground state of 
patterned magnetic dots, the dynamics of vortices have 
attracted considerable attention. Being displaced from 
its equilibrium position in the dot center, the vortex 
core reveals sub-GHz frequency oscillations with a nar- 
row linewidth [2,7,12]. The oscillations of the vortex 
core are governed by a competition of the gyroforce, 
Gilbert damping force, spin transfer torque, and restor- 
ing force. The restoring force is determined by the vor- 
tex confinement in a nanodot. Vortex core oscillations 
with small amplitude can be well described in the linear 
regime, but for increasing of the STNO output power, a 
large-amplitude motion has to be excited. In the regime 
of large-amplitude spin transfer-induced vortex gyr- 
ation, it is important to take into account nonlinear 
contributions to all the forces acting on the moving vor- 
tex. The analytical description and micromagnetic simu- 
lations of the magnetic field and spin transfer-induced 
vortex dynamics in the nonlinear regime have been pro- 
posed by several groups [12-22], but the results are still 
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Figure 1 Magnetic vortex dynamics in a thin circular FeNi nanodot. Vortex core steady-state orbit radius UqW in the circular FeNi nanodot of 
thicl<ness i = 7 nm and radius fi= 100 nm vs. current J perpendicular to tlie dot plane. Solid black lines are calculations by Equation 7; red circles 
mark the simulated points. Inset: sketch of the cylindrical vortex state dot with the core position X and used system of coordinates. 



contradictory. It is unclear to what extent a standard 
nonlinear oscillator model [13] is applicable to the vor- 
tex STNO, how to calculate the nonlinear parameters, 
and how the parameters depend on the nanodot sizes. 

In this paper, we show that a generalized Thiele approach 
[23] is adequate to describe the magnetic vortex motion in 
the nonlinear regime and calculate the nanosize vortex core 
transient and steady orbit dynamics in circular nanodots 
excited by spin-polarized current via spin angular momen- 
tum transfer effect. 

Methods 

Analytical method 

We apply the Landau-Lifshitz-Gilbert (LLG) equation of 
motion of the free layer magnetization ih = -ym x Heff + 
ctGin X rh + yXs , where m = M/M^, is the saturation 
magnetization, > 0 is the gyromagnetic ratio, Hgff is the 
effective field, and is the Gilbert damping. We use a 
spin angular momentum transfer torque in the form 
suggested by Slonczewski [24], = o/m x (m x P), 
where a= ^/{2\e\LMs), tj is the current spin polarization 
[t] = 0.2 for FeNi), e is the electron charge, P is direction of 
the reference layer magnetization, and / is the dc current 
density. The current is flowing perpendicularly to the 
layers of nanopillar and we assume P = Pz. The free layer 
(dot) radius is R and thickness is L. 

We apply Thiele's approach [23] for the magnetic vortex 
motion in circular nanodot (inset of Figure 1). We assume 
that magnetization distribution can be characterized by 
a position of its center X(i) that can vary with time 
and, therefore, the magnetization as a function of the 



coordinates r and X(t) can be written as m(r,t) = m(r,X(i)). 
Then, we can rewrite the LLG equation as a generalized 
Thiele equation for X(f): 



(1) 



where W is the total magnetic energy, a,/S = x,y, and 
da = (^/t^Xa- The components of the gyrotensor G, damping 
tensor D, and the spin-torque force can be expressed as 
follows [16]: 



Ga^(X) =— J d^r(dam x d^m) ■: 
Dafi{X) = -Ug — J d^r[dam x dp 
f^j,(X) = MsCrJP-J <i^r(m x a^m; 



m 



(2) 



We assume that the dot is thin enough and m does 
not depend on z-coordinate. The magnetization m{x,y) 



has the components mx 



im. 



■y ' 



2w/{l- 



WW 



and 



= {l-ww)/{l + ww) expressed via a complex function 
^(CC) [25]. Inside the vortex core, the vortex configur- 
ation is described as a topological soliton, w(C, C) = /(C)' 
[/(O I ^ 1, where fX) is an analytic function. Outside the 
vortex core region, the magnetization distribution is w 
=/(()/[/-(f)|, [/(Ol >1. For describing the vortex 
dynamics, we use two-vortex ansatz (TVA, no side surface 
charges induced in the course of motion) with function /(() 
being written as /(C) = -«C(C-5)(5C-l)/c(l + \s\^) [26], 
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where C is the vortex chirahty, (= {x + iy)IR, s = Sx + isy, 
s = X/7?, c = RJR, and Rc is the vortex core radius. 

The total micromagnetic energy W = W^^ + 
Wex + in Equation 1 including volume and sur- 
face magnetostatic energy, exchange W^.^ energy, 
and Zeeman Wz energy of the nanodot with a displaced 
magnetic vortex is a functional of magnetization distri- 
bution \F[m(r, t)]. Using m = in(r, X(f)) and integrating 
over-the-dot volume and surface, the energy W can be 
expressed as a function of X within TVA [16]. The Zeeman 

energy is related to Oersted field H; = (0,///,0) of the 

spin-polarized current, WzOQ = - dVm{r, X) • H/. We 
introduce a time-dependent vortex orbit radius and phase 
by s = M exp(i<t>). The gyroforce in Equation 1 is determined 
by the gyrovector G = Gz, where G=Gz = G^y The func- 
tions G(s) and W{s) depend only on m = |5| due to a circular 
symmetry of the dot. G(0) = InpM^L/y, where p is the vor- 
tex core polarity. The damping force DX and spin-torque 
force FsT-are functions not only on m = |5| but also on direc- 
tion of s. Nonlinear Equation 1 can be written for the circu- 
lar dot in oscillator-like form 

is + (x)g{u)s = -d{u)s + ix{u)s-id„{s)s, (3) 

where coeiu) = {R^u\G{u)\y ^dW{u)/du is the nonlinear 
gyrotropic frequency, d{u) = - D{u)/\G(u)\ is the nonlinear 
diagonal damping, D = D^ = Dyy, dn{s) = - D^{s)/\G{u)\ is 
the nonlinear nondiagonal damping, and/(M) = a{u)/\G{u)\. 
It is assumed here that FsAs) = a(u){z x s) [14], where 
a is proportional to the CPP current density / and a 
(0) = jiRLM.aJ. 

To solve Equation 3, we need to answer the following 
questions: (1) can we decompose the functions W{s), G 
(s), Dafiis), and Fst^s) in the power series of m = \s\ and 
keep only several low-power terms? and (2) what is the 
accuracy of such truncated series accounting that m = |5| 
can reach values of 0.5 to 0.6 for a typical vortex STNO? 
Some of these functions may be nonanalytical functions 
of M = |5|. If the answer to the first question is yes, then 
we should decompose W{s) up to m*, Fst^s) up to u^, 
and G(s), Da^is) up to w^-terms to get a cubical equation 
of the vortex motion. The series decomposition of G(s) 
does not contain w^-term; it contains only small c^M^-term, 
G(m) = G(0)[1 - O(c^M^)], although G(m) essentially decreases 
at large u, when the vortex core is close to be expelled firom 
the dot [16]. The result of power decomposition of the total 
energy density w(m) = W{u)/M'^V is 

w{u) = w{0) +l^KU^ + ^k'u'^, M=|5|, (4) 
2 4 

and the coefficients are 



0 

lK{fi,RJ) = 2nj d/-^[ll{t)-I{t)h{t)] +\(^-^' 

0 

n ICR 

^ 15c Ms ' 

where /(f) = / dppj,{tp), h{t) = f dpp-^l-p'fj.itp), 
Jo Jo 

h{t)= j dp{l + py2{tp), l^ = L/R, Le = V2A/Ms, and 
Jo 

q = l + l5{[n2-l/2)RJ8R. 

There is an additional contribution to k/2, 2{LJRf, 
due to the face magnetic charges essential for the nano- 
dots with small R [27]. The contribution is positive and 
can be accounted by calculating dependence of the equilib- 
rium vortex core radius (c) on the vortex displacement. 
This dependence with high accuracy at cm << 1 can be 
described by the function c(m) = c(0)(l - m^)/(1 + m^). 
Here, c(0) is the equilibrium vortex core radius at s = 0, 
for instance c(0) = 0.12 {Rc = 12 nm) for the nanodot 
thickness L = 7 nm. 

The nonlinear vortex gyrotropic frequency can be 
written accounting Equation 4 as 

toeiu) = coo[l+Nu^], (5) 

where the linear gyrotropic frequency is <ao = jMsKifi, R,J)I2, 
and N(fi, R) = xifi, R)lK{fi, R). 

The frequency Wq = yMsK{f^,R, 0)/2 was calculated in 
[26] and was experimentally and numerically confirmed 
in many papers. The nonlinear coefficient N(fi,R) depends 
strongly on the parameters yS and R, decreasing with ^ and 
R increasing. The typical values of N{li,R) at / = 0 are equal 
to 0.3 to 1. 

The last term in Equation 3 prevents its reducing to a 
nonlinear oscillator equation similar to the one used for 
the description of saturated STNO in [13]. Calculation 
within TVA yields the decomposition dn{s) = <i° + d\sxSy, 
where (i° = 0 , i.e., the term containing dn{s) ~ UqiP' <<1 
can be neglected. Then, substituting s = u exp(iO) to 
Equation 3, we get the system of coupled equations 

^-Ug{u) = d{u) — , M = [j(M)-(^(M)<i>] M. (6) 

Equation 3 and the system (6) are different from the 
system of equations of the nonlinear oscillator approach 
[13]. Equations 6 are reduced to the autonomous oscilla- 
tor equations u / u = x{u)-d{u)o}Q{u) and <i> = (jOg{u) 
only if the conditions d^ < < \ and dx< < are satisfied 
and we define the positive/negative damping parameters 
[13] as r^{u) = d{u)(OQ{u) and r_{u) = j(m). We note that 
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reducing the Thiele equation (1) to a nonlinear oscillator 
equation [13] is possible only for axially symmetric 
nanodot, when the functions W{s), G{s), d{s) and xis) 
depend only on m= |s| and the additional conditions 
d„ < < 1, < < 1, and dx< < Wg are satisfied. The non- 
linear oscillator model [13] cannot be applied for other 
nanodot (free layer) shapes, i.e., elliptical, square, etc., 
whereas the generalized Thiele equation (1) has no 
such restrictions. 

The system (6) at m = 0 yields the steady vortex oscil- 
lation solution UoU) > 0 as root of the equation j(mo) = d 
(mo)wg(wo) for xiO) > d{0)cjo (J > /^i) and Mq = 0 otherwise. 
If we use the power decompositions cociu) = wo + cjiu^, 
d{s) = do + diu^ , and j(m) = Jo + JiM^ for the nonlinear 
vortex frequency, damping, and spin-torque terms, re- 
spectively, and account that the linear vortex frequency 
contains a contribution proportional to the current density 
Wo(/) = + ^el> where (Ug = {8n/15){YR/c)q [12,16], then 
we get the vortex core steady orbit radius at / > 

"o(/) =A(/)v'///ci-l, A'a) 

[diMD+doMD-XiiJ)]' 

The model parameters are Jd = docj^/ {ya /2-doCJe) , 
do = aG[5 + 4 ln{R/Rc)]/8, d^ = llaa/e, xo = ya-J/2. The ratio 
Xi^Xo = 0{c'i/) < < 1, therefore, the nonlinear parameter 
can be neglected. The statement about linearity of the 
ST-force agrees also with our simulations and the 
micromagnetic simulations performed in [12,19]. The 
coefficient X(J) describes nonlinearity of the system 
and decreases smoothly with the current / increasing. 

Numerical method 

We have simulated the vortex motion in a single 
permalloy (Fe2oNi8o alloy, Py) circular nanodot under 
the influence of a spin-polarized dc current flowing 
through it. Micromagnetic simulations of the spin-torque- 
induced magnetization dynamics in this system were 
carried out with the micromagnetic simulation package 
MicroMagus (General Numerics Research Lab, Jena, 
Germany) [28]. This package solves numerically the LLG 
equation of the magnetization motion using the optimized 
version of the adaptive (i.e., with the time step control) 
Runge-Kutta method. Thermal fluctuations have been 
neglected in our modeling, so that the simulated dynamics 
corresponds to T= 0. Material parameters for Py are as 
follows: exchange stiffness constant A = 10"^ erg/cm, 
saturation magnetization Ms = 800 G, and the damping 
constant used in the LLG equation ag = 0.01. Permalloy 
dot with the radius R = 100 nm and thickness L = 5,7, and 
10 nm was discretized in-plane into 100 x 100 cells. No 
additional discretization was performed in the direction 
perpendicular to the dot plane, so that the discretization 



cell size was 2 x 2 x L nm''. In order to obtain the vortex 
core with a desired polarity (spin polarization direction of dc 
current and vortex core polarity should have opposite direc- 
tions in order to ensure the steady-state vortex precession) 
and to displace the vortex core from its equilibrium position 
in the nanodot center, we have initially applied a short 
magnetic field pulse with the out-of-plane projection 
of 200 Oe, the in-plane projection = 10 Oe, and the 
duration At = 3 ns. Simulations were carried out for 
the physical time t = 200 to 3,000 ns depending on the 
applied dc current because for currents close to the 
threshold current /d, the time for establishing the vortex 
steady-state precession regime was much larger than for 
higher currents (see Equation 8 below). 

Results and discussion 

Calculated analytically, the vortex core steady orbit radius 
in circular dot Uoif) as a function of current / is compared 
with the simulations (see Figure 1). There is no fitting ex- 
cept only taking the critical current /d value from simula- 
tions. Agreement is quite good, confirming that all the 
nonlinear parameters of the vortex motion were accounted 
correctly. The steady orbit radius Mo(/) allows finding 
the STNO generation frequency cjg(/) = Wo(/) + (jOiU^{)), 
which increases approximately linearly with / increasing up 
to the second critical current value Jc2 when the steady 
oscillation state becomes unstable (see Figure 2). The in- 
stability is related with the vortex core polarity reversal 
reaching a core critical velocity or the vortex core expelling 
from the dot increasing the current density / [12,16]. We 
simulated the values of Jc2 = 2.7, 5.0, and 10.2 MA/cm for 
the dot thickness L = 5, 7, and 10 nm, respectively. The 
calculated STNO frequency is 15 to 20% higher than 
the simulated one due to overestimation of Wq within 
TVA for yS =0.1. The calculated nonlinear frequency part 
is very close to the simulated one, except the vicinity of 
/c2, where the analytical model fails. 

Our comparison of the calculated dependences MqC/) 
and cjgU) with simulations is principally different from 
the comparison conducted in a paper [19], where the au- 
thors compared Equations 5 and 7 with their simulations 
fitting the model-dependent nonlinear coefficients N and 
A from the same simulations. One can compare Figures 1 
and 2 with the results by Grimaldi et al. [20], where the 
authors had no success in explaining their experimental 
dependences Uo(J) and og(J) by a reasonable model. The 
realistic theoretical nonlinear frequency parameter N for 
Py dots with L = 5 nm and R = 250 nm should be larger 
than 0.11 that the authors of [21] used. TV =0.25 can be 
calculated from pure magnetostatic energy in the limit 
yS — > 0 (inset of Figure 2). Accounting all the energy contri- 
butions in Equation 4 yields N = 0.36, which is closer to 
the fitted experimental value AT = 0.50. 
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Current density J (MA/cm^) 

Figure 2 The vortex steady-state oscillation frequency vs. current. The nanodot thickness /_ is 5 nm (1), 7 nm (2), and 10 nm (3), and radius 
is fi = 1 00 nm. The frequency is shown within the current range of the stable vortex steady-state orbit, J^, <J< J^2- Solid black lines are calculations 
by Equation 5; red squares mark the simulated points. Inset: the nonlinear vortex frequency coefficient vs. the dot thickness for ft = 100 nm 
and J = 0 accounting all energy contributions (1) and only magnetostatic contribution (2). 



The system (6) can be solved analytically in nonlinear 
case. Its solution describing transient vortex dynamics is 



"o(/) 



1(2(0) " 



1 exp 



(8) 



where m(0) is the initial vortex core displacement and 
l/r+(/) = 2dQ(0Q{J/Jci-l) is the inverse relaxation time 
for />/ci (order of 100 ns). u{tjci)'^l/yft at t ^ oo and 
/ = Jci- If / < /ci, the orbit radius u{t, ]) decreases expo- 
nentially to 0 with the relaxation time 1/t-{J) = cIqUq 
|///ci~l|- The divergence of the relaxation times r+ at 
/ = Jci allows considering a breaking symmetry second- 
order phase transition from the equilibrium value Mq = 0 
to finite Uo{})'^\/j /}c\-^ defined by Equation 7. Equa- 
tions 7 and 8 represent mean-field approximation to the 
problem and are valid not too very close to the value of /= 
Jci, where thermal fluctuations are important [13,21]. 

Equation 8 describes approaching of the vortex orbit 
radius to a steady value Mo(/) under influence of dc spin- 
polarized current. We note that the corresponding relax- 
ation time is determined by only linear parameters, whereas 
the orbit radius (7) depends on ratio of the nonlinear 
and linear model parameters. The solution of Equation 8 
is plotted in Figure 3 as a function of time along with 
micromagnetic simulations for circular Py dot with 



thickness L = 7 nm and radius R = 100 nm. The vortex 
was excited by in-plane field pulse during approximately 
the first 5 ns, and then the vortex core approached the sta- 
tionary orbit of radius MoO- We estimated m(0) after the 
pulse as m(0) = 0.1 and plotted the solid lines without any 
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Figure 3 Instant vortex core orbit radius vs. time for different 
currents. The results are within the current range of the stable 
vortex steady-state orbit, J^, <J<J^2 (5.0 MA/cm^). The nanodot 
thickness is /_ = 7 nm and the radius is R= 100 nm. Solid lines are 
calculations of the vortex transient dynamics by Equation 8, and 
symbols (black squares, red circles, green triangles, and blue rhombi) 
mark the simulated points. 
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fitting except using the simulated value of the critical 
current fd- Overall agreement of the calculations by 
Equation 8 and simulations is quite good, especially for 
large times t > 3t+, although the calculated relaxation time 
T+ is smaller than the simulated one due to overestimation 
of cOq within TVA. The typical simulated ratio /c2/fci ~ 1-5; 
therefore, minimal r+ ~ 20 to 30 ns. But the transient time 
of saturation of u{t,J) is about of 100 ns and can reach 
several microseconds at J/Jd < 1.1. The simulated value of 
A = 0.83, whereas the analytic theory based on TVA yields 
the close value of = 0.81. 

Typical experiments on the vortex excitations in nanopil- 
lars are conducted at room temperature T= 300 K without 
initial field pulse, i.e., a thermal level m(0) should be suffi- 
cient to start vortex core motion to a steady orbit. To find 
the thermal amplitude of m(0), we use the well-known 
relation between static susceptibility of the system x and 
magnetization fluctuations (Af^)y-(/Wa)^= (/f^r/l/)/^^. 
The in-plane components are Xxx ~ Xyy = {^^/'^) ^ 
{yMs/cj'o), and M = ^M,s, where ^=2/3 within TVA [26]. 
This leads to the simple relation (u^)j- = [ksT/M^V) 
(yMs/ cJq) ■ It is reasonable to use Ut{0) = \/ {u'^)j for in- 
terpretation of the experiments. Mr(0) ~ 0.05 (5 nm in ab- 
solute units) for the dot made of permalloy with L = 7 nm 
andJ? = 100nm. 

The nonlinear frequency coefficient N(fi,R,J) = K'(Ji,R,J)/ 
K(fi,R,J) is positive (because of k,k' >0 for typical dot 
parameters), and it is a strong function of the dot geomet- 
rical sizes L and R and a weak function of /. For the dot 
radii R>>L^, N(fi, R, 0) « 0.21 - 0.25 (the magnetostatic 
limit, see inset of Figure 2). If 7? > > L^, and /S — > 0, then N 
iji, R, 0) ~ 0.25 [14]. For the realistic sizes of free layer in a 
nanopUlar [R is about 100 nm and L = 3 to 10 nm), this 
coefficient is essentially larger due to finite ji and exchange 
contribution, and it can be of order of 1. The exchange 
nonlinear contribution /c'gx is important for R < 300 nm. 
However, the authors of [19-21] did not consider it at all. 
Note that M0.089, 300 nm, 0) « 0.5 recently measured [29] 
is two times larger than 0.25. The authors of [19] suggested 
to use an additional term ~u in the magnetic energy fitting 
the nonlinear frequency due to accounting a u -contribu- 
tion (Af=0.26) that is too small based on [14], while the 
nonlinear coefficient N(p,R) calculated by Equation 5 for 
the parameters of Py dots {L = 4.8 nm, R = 275 nm) [19] is 
equal to 0.38. Moreover, the authors of [19] did not account 
that, for a high value of the vortex amplitude u = 0.6 to 0.7, 
the contribution of nonlinear gyrovector G{u) °< c^u^ to the 
vortex frequency is more important than the u -magnetic 
energy term. The gyrovector G{u) decreases essentially for 
such a large u resulting in the nonlinear frequency increase. 
The TVA calculations based on Equation 5 lead to the 
small nonlinear Oe energy contribution /c'oe> whereas 



Dussaux et al. [19] stated that Ar'oe is more important 
than the magnetostatic nonlinear contribution. 

Conclusions 

We demonstrated that the generalized Thiele equation 
of motion (1) with the nonlinear coefficients (2) considered 
beyond the rigid vortex approximation can be successfully 
used for quantitative description of the nonlinear vortex 
STNO dynamics excited by spin-polarized current in a 
circular nanodot. We calculated the nonlinear parameters 
governing the vortex core large-amplitude oscillations and 
showed that the analytical two-vortex model can predict 
the parameters, which are in good agreement with the ones 
simulated numerically. The Thiele approach and the energy 
dissipation approach [12,19] are equivalent because they 
are grounded on the same LLG equation of magnetization 
motion. The limits of applicability of the nonlinear oscilla- 
tor approach developed for saturated nanodots [13] to vor- 
tex STNO dynamics are established. The calculated and 
simulated dependences of the vortex core orbit radius u{t) 
and phase 0(f) can be used as a starting point to consider 
the transient dynamics of synchronization of two coupled 
vortex ST nano-oscillators in laterally located circular 
nanopillars [30] or square nanodots with circular nano- 
contacts [31] calculated recently. 
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